Critical properties of a three dimensional p-spin model 
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In this paper we study the critical properties of a finite 
dimensional generalization of the p-spin model. We find evi- 
dence that in dimension three, contrary to its mean field limit, 
the glass transition is associated to a diverging susceptibility 
(and correlation length). 



I. INTRODUCTION 

Two different transition mechanisms are known in spin 
glass mean-field theory |l],^], according to the form of 
the random Hamiltonian. In some models, like the infi- 
nite range Shcrrington-Kirkpatrick (S-K) model for spin 
glasses, there is a second order glassy transition with di- 
verging spin-glass susceptibility and continuous replica 
symmetry breaking. In other models, whose prototype is 
the Random Energy Model, the transition is first order 
with a Gibbs-Dimarzio like entropy crisis. Other exam- 
ples of models with the second type of behaviour are spin 
models with p-spin interaction, with p > 2, both for Ising 
spins and for spherical spins. 

Numerical simulations indicate that the first type of 
mechanism describes the ergodicity breaking transition of 
finite dimensional spin glasses || . The second mechanism 
is more appropriate to describe the behavior of structural 
glasses j|. The passage from mean-field to finite dimen- 
sion is in both cases highly non trivial. Despite many 
progresses || the problem of including fluctuations in the 
description of the finite dimensional spin glass transition 
is far from being achieved. For that reason the test of 
the mean-field picture has been left in the last 15 years 
to the numerical study of the Edwards- Anderson model, 
which admits the SK as infinite dimensional limit. 

Strangely enough there are only few numerical studies 
of finite dimensional spin models that could have a tran- 
sition homologous to the mean-field discontinuous transi- 
tion. Given the relevance of this transition to structural 
glasses, the study of such finite dimensional models is of 
primary importance. 

Up recentlyto our knowledge, the only studies ap- 
peared in the literature, are these references || and 0. 
In reference 0] it was proposed a generalization of the p- 
spin model which presents a phenomenology reminiscent 
to that of structural glasses. However, the difficulty to 



decide about the existence of a phase transition and the 
presence of a spurious symmetry, makes necessary to re- 
sort to better conceived models without spurious effects. 

In this paper we introduce, and study numerically, a 
finite dimensional model with ./V spin per sites and p-spin 
interactions, that for all dimensionality tends to a mean- 
field behaviour for N — > oo. As it happens in ordinary 
field-theory, this large N approach is complementary to 
the high-dimensionality approach. We study this model 
for p = 4 in D = 3. A complementary study of the same 
model (still for p — 4 and D — 3) in the low temperature 
regime can be found in ref. Some numerical simula- 
tion of the p = 3 model for D — 4 can be found in ref. 



II. A FINITE DIMENSIONAL P-SPIN MODEL 



The long range p-spin model 
Hamiltonian is given by 
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11 Ol is soluble and its 



(1) 



where the variables J are random with zero average 

and variance \/N { p-^ (the same result is obtained for 

Gaussian distributed variables and in the case J = 
±l/Ar((p-i)/2)) and the 

spins are Ising variables (also the 
spherical case is soluble). 

This Hamiltonian can be generalized in many ways in 
finite D. The way we follow in this paper is the follow- 
ing: we consider a D dimensional square lattice with N 
Ising spins S" (a = 1, N) in each site x of the lattice. 
For any given couples of nearest neighbour sites there are 
(2N)\ / p\(2N — p)\ possible groups g of p-spins. We con- 
sider the product of all the spins in each group, and we 
couple them with a random variable J g . The resulting 
Hamiltonian, with transparent notation, is 



H 



e e ^n^ 



(2) 



<x,y> g e(x,y) 



The J p are chosen 



where we have relabeled the spins 
independently from link to link and group to group and 
are equal to ±1 with probability 1/2. The mean field 
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limit is recovered both for high dimension (D — > oo) and 
finite N or for large N and finite D. Indeed it is possi- 
ble to construct a loop expansion for the development in 
powers of 1/N jnj. In this paper we present a numerical 
study of the model in D = 3 for low values of N above 
the transition point. 

We have simulated the physics of the model through 
Monte Carlo method with the Metropolis algorithm. We 
have studied in detail the cases p — 4 and N = 3 and 
4. As we will see the results do not seem to depend 
qualitatively on N in this range of N (for N = 2 and 
p = 4 the model is isomorphic to the usual short range 
Edwards Anderson model for spin glasses). 

In this first study of the model we discuss mainly three 
issues: 

• The existence of a glass transition by means of the 
study of the thermodynamics of the model, through 
the behaviour of the energy and the entropy in sim- 
ulations of "cooling experiments" . 

• The behaviour of the time dependent auto- 
correlation function at equilibrium and the grow- 
ing of the relaxation time as the glass transition is 
approached. 

• The existence of a growing spin-glass correlation 
length. 

Let us start with the discussion of cooling experiments. 
In figure [l] and ^| we show the energy as a function of the 
temperature for different cooling rates for N = 3 and 4 
respectively. The cooling rate k is equal to the inverse 
of the number of Monte Carlo sweeps done at each tem- 
perature. We recognize in both cases the typical curves 
of systems undergoing a glass transition, and remaining 
frozen below a cooling rate dependent freezing temper- 
ature. In the figure with N — 4 we have plotted for 
comparison purposes the line corresponding to the first 
term of the high temperature expansion. We see that 
until quite near to the freezing the energy of the system 
remains close to that line. This is reminiscent to what 
happens in mean-field, where the line is followed up to 
the transition point. 

In figure [| we study the dependence of the energy on 
the cooling rate for the N = 4 model for T — 2 and 
T = 4 as a function of the inverse cooling rate. The 
data are compatible with a power law relaxation of the 
kind E(k) = E^ + Ak u with a temperature dependent 
exponent u. For instance a fit of the data for T = 2, 4 
gives E(k)\ t =2 = -65.8 + 11.6 x k 021 and E(k)\ t =4 = 
—65.1 + 11.9 x K - 23 . This dependence contrasts with 
the much slower logarithmic dependence observed in real 
glasses and is the first sign of criticality in the system. 

In the equilibrium regime (which is reached for suf- 
ficient long simulations) we can obtain free-energy and 
entropy integrating the data of the internal energy and 



taking into account that at infinite temperature the en- 
tropy per spin is S(T — oo) = log 2. The free-energy is 
reconstructed as 

F = log2-T [ f E(J3)d0, (3) 
Jo 

and from it the entropy, that we plot in figure [|. We ob- 
serve that S(T) behaves linearly in a wide range of tem- 
peratures suggesting the validity of the Gibbs-Di Marzio 
transition mechanism for this model. As we will see the 
model with = 3 seems to have a transition around 
T = 2.6 while the linear extrapolation of the entropy 
vanishes only at T w 1.8. However it is not the total 
entropy, but the "configurational entropy" (associated to 
the number of possible metastable states) that should 
vanish at the transition. 

We now turn to the more difficult question of the iden- 
tification and characterization of the phase transition in 
the model. We have studied that issue limiting ourselves 
to the case = 3. The quantity over which we have 
concentrated is the overlap correlation length. We sim- 
ulated two identical replicas in parallel (a and r), and 
after thermalization we measured the overlap correlation 
function: 

V N N 

G ^ = jpv E E E <°r +X T?T? +x . (4) 

The spin glass susceptibility is defined as 

Xsg = J d 3 xG(x). (5) 

The data for the function G(x) are shown in figure |5| 
for T = 3, together with the the best fit of the form 
G(x) = A/(x + l) 1+r >exp(-x/£) with r? = 0. We have 
done similar fits at different temperatures and in this way 
we have extracted the value of the correlation length. 
If we plot this correlation length as a function of the 
temperature (fig. |J) we see that the data are best fitted 
by the power form f « (T - T c ) v ', with T c = 2.62 and 
v = 0.71. However from figure we see that the data are 
also well compatible with v = 2/3 as one could expect 
from scaling arguments (see below). 

Differently from mean field there is a correlation length 
growing in the system, that suggests a second order phase 
transition. 

Similar conclusions also come from the analysis of the 
spin glass susceptibility, which grows of about a factor 
ten in the range 3 < T < 4.4. The data, shown in figure 
U are roughly compatible with a power law divergence of 
the susceptibility as ~ £ 2 4 , which using the value 2/3 for 
v corresponds to \sg ~ (T — T c ) 7 with 7 = 1.6. These 
values for the critical exponents 7 and v are definitely 
different from those of the Edwards Anderson spin glass 
models: they are a factor 2-3 times smaller |1. 
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The last aspect that we have studied is the equilib- 
rium relaxation, and the relation of the relaxation time 
with the correlation length. In figure || we show the time 
autocorrelation function 

c W = ^I>, ? (*K ? (0) (6) 

for various temperature for N = 4. We see that a form 
C(t) = Acxp(— (t/r) 13 ) fits excellently the data. We ex- 
tract from that the relaxation time r (figure |l^) and the 
exponent [3 that we plot in figure [ill 

We also done the same analysis for N = 3 with similar 
results. It is interesting to plot r versus £ which shows 
the compatibility of our data with the relation r ~ £ 2 
(see figure [jj]). This scaling form is at variance with 
the results of JtJ for the other short range p-spin model 
we mentioned in the introduction, and experiments in 
structural glasses. The high value of z we find. i.e. z = 8 
(which by coincidence is quite similar to the value in the 
Edwards Anderson spin glass model) implies a violent 
divergence of the correlation time near the transition and 
it is responsible of the large correlation time needed to 
reach thermalization. 

The difference of behaviour with respect to the mean- 
field can be rationalized with the argument presented in 
the next section, which also implies v = 2/3. 

III. A POSSIBLE THEORETICAL 
INTERPRETATION 

The difference among the short range model and the 
homologous infinite range model (which can be solved in 
the mean field approximation) are quite striking. No pre- 
cursor signs of the transition are present in the infinite 
range model and the spin glass susceptibility remains fi- 
nite up to the transition point. Here we will present and 
argument which suggests that in short range models with 
quenched local disorder the spin glass susceptibility is di- 
vergent as in second order phase transitions. 

The precise reasons for this discrepancy are not clear 
to us. The simplest scenario we have considered is the 
following. The disorder induces fluctuations in the lo- 
cal transition temperature. For a given region of radius 
R centered around the point x we can define and effec- 
tive critical temperature Tr(x). It is natural to assume 
that the x-dependent fluctuation of Tr(ie) is a quantity 
of order R~ D / 2 , i.e. 

T R {x) = T c + 5T R (x) (7) 

with (STr(x) 2 ) oc R D . At a given temperature T the 
regions of size R such that Tr(x) > T are strongly cor- 
related. The typical radius of these region increases as 
(T - T c fl D suggesting therefore that v = 2/D. This 
value of v is peculiar for random systems. Indeed there 



are general arguments that show that for second order 
phase transitions disorder is relevant if v > 2/D [[l2| . 

We also notice that the same value of v can be ob- 
tained if we assume that the specific heat has a discon- 
tinuity at the phase transition as predicted by the mean 
field analysis. In other words we suppose that the specific 
heat exponent a is equal to zero. The usual scaling law 
a = 2 — Dv implies the result v = 2/D. This argument 
is suggestive, but the coincidence its prediction with the 
value of v we find may be fortuitous of the dimension 
three. An investigation of the model in higher dimen- 
sions will give some information on the validity on this 
conjecture (numerical simulations in four dimensions for 
the p = 3 model Q suggest that in this case v is around 
.5). 

IV. CONCLUSIONS 

In this paper we have studied by Monte Carlo a fi- 
nite dimensional version of the p-spin model. We find a 
scenario for freezing that mixes typical features of struc- 
tural glasses, like strong cooling rate dependence of the 
low temperature energy, with features of second order 
phase transitions with power law growing of the corre- 
lation length and critical dynamics. This is a genuine 
finite dimensional effect due to the quenched disorder 
which can be understood qualitatively with the argument 
we have given. A full theoretical understanding should 
come from the inclusion of non-perturbative effects in the 
theory. The simulations of disordered finite dimensional 
analogous of systems with "one step replica breaking" is 
just at the beginning and much progress can be expected 
in the future. 
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FIG. 1. The energy for N = 3 and three different values 
of the cooling rate k = 2~ 5 , 2" 8 , 2" 17 . 




FIG. 2. The energy for N = 4 and three different values of 
the cooling rate (2~ 5 , 2~ 7 , 2~ 12 ), together with the first term 
of the high temperature expansion. The curves stay close to 
the first term of the high temperature expansion until very 
close to where they fall off equilibrium. 




FIG. 3. The energy for N — 4 as a function of the inverse 
cooling rate. From top to bottom T = 8, 6, 4, 2. The full lines 
are the power law fits discussed in the text. 
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FIG. 4. The entropy as a function of the temperature for 
N = 3. 
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FIG. 5. Overlap correlation function as a function of dis- 

r »r o m o i -,i ,i c , ,., r FIG. 7. Correlation length t to the power -2/3 versus 

tance for TV = 3, 1 — 3, together with the ht or the form ,. „ 



c ( :r ) = ex P(— x /0- The value of £ from the fit is £ = 2.5. 



temperature for N = 3. The curve is a linear fit. 
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FIG. 6. Inverse of the correlation length versus tempera- 
ture for TV = 3 and ht £ _1 = (T - T c ) v , v = 0.71, T c = 2.62. 



FIG. 8. Spin glass susceptibility xsg versus £ for TV = 3 
and the power best fit which gives x ~ £ 2 ' 4 
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FIG. 9. Time dependent auto-correlation function at 
equilibrium in the N — 4 model for, from top to bottom, 
T = 5, 6, 7, 8, 9, 10 the lines are stretched exponential fits. 



FIG. 11. Stretching exponent versus temperature for 
N = 4. The line is a linear fit of the region T < 10: 
/3 = T x 0.093 - 0.29. 




FIG. 10. Relaxation time extracted from the autocorre- 
lation function as a function of the temperature for N = 4. 
The line is a power law fit t = 117.8 X (T - 6.4) 2 .1. 



FIG. 12. Relaxation time versus £ for N = 3. The line is 
a power law fit with power z = 8. 



